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Abstract 

Primary graviton spectra, produced via stimulated emission from an initial Bose-Einstein 
distribution, are enhanced for typical scales larger than the redshifted thermal wavelength. 
A mixed state of phonons induces a secondary graviton spectrum which is hereunder com- 
puted in terms of three parameters (i.e. the number of phonon species, the tensor-to-scalar 
ratio and the thermal wavelengths of the mixture). The primary and secondary graviton 
spectra are shown to be sensitive, respectively, to the first-order and second-order correlation 
properties of the initial quantum mixture so that the semiclassical theory is argued to be 
generally inadequate in this context. For particular values of the parameters the secondary 
contribution may turn out to be comparable with the primary spectrum over large-scales. 



Electronic address: niassinio.giovannini@cern.ch 



1 Introduction 



The analysis of the Cosmic Microwave Background (CMB) anisotropies and polarization 
[H El El H] , of the extended galaxy surveys [5l E] and of the supernova data sets [71 18] aims 
at a highly accurate statistical determination of the cosmological parameters belonging to 
a class of simplified scenarios developed around the so-called concordance paradigm also 
dubbed, sometimes, ACDM modefl In the concordance paradigm the tensor modes are 
assumed to be absent but there are good reasons to believe in their existence. It is plausible 
to think that, in the near future, we shall be forced to deal with some sort of TACDM 
scenario [S] (with T standing for the tensor modes) where experiments operating at high- 
frequencieqf|(such as the wide-band interferometers [IHIITT] or devices operating at even larger 
frequencies [121 US]) and CMB polarization experiments could be jointly analyzed to probe 
possible excesses both over large length-scales (i.e. small frequencies) and over comparatively 
higher frequencies. 

The precise nature of the pre-inflationary initial conditions calls for supplementary as- 
sumptions [13] and it is anyway not explained in the framework of the simplistic ACDM 
paradigm. It is useful to speculate that the spectrum of relic gravitons has a steeper slope 
at large scales. This is what happens if thermal corrections to the graviton spectrum are 
consistently included in the scenario. The thermal corrections arise if a Bose-Einstein dis- 
tribution of spin-two fluctuations of the geometry was present around 65 efolds prior to the 
end of a conventional inflationary phase (see, e.g. [151 HI HI])- The correction to the power 
spectrum takes a standard form which can be written as [TH HH [201 1211 |22l |23l |23] 

P«(A;,M = coth(A^PT(fc), (1.1) 

where = T is the comoving thermal wavenumber in units h = c = kb = ^] -Pt(^) is the 
spectrum in the absence of thermal corrections so that the two spectra are not distinguishable 
when k ^ k'Y while the spectral slopes differ for wavelenghts A > At = 27r/A;x- In the limit 
A ^ At, Pt^''(A;, ^t) — (^/^t)^^ -Pt(^)- More general situations can be contemplated if the 
initial state of the gravitons is mixed but not exactly thermal. In all these situations particles 
are produced by stimulated emission as opposed to the spontaneous emission operating when 
the initial state is the vacuum. For wavenumbers k S> k^, as a consequence of Eq. (11.11) . the 
stimulated and spontaneous emissions lead, in practice, to the same results. 

^In the acronym the A term parametrizes the dark energy component and CDM stands for the cold dark 
matter component 

■^Current CMB and large-scale observations are commonly phrased in terms of wavenumbers expressed 
in units of Mpc~^ while the operational windows of wide-band interferometers (cross-correlated among 
themselves and with some bar detector) are given in Hz. It is useful to recall, in this respect, that the pivot 
scale at which the CMB power spectra are assigned corresponds to kp = 0.002 Mpc^^ (related to an effective 
multipole £p — 30). In the natural units used throughout this paper (i.e. h = c = kb = 1) the frequency 
lyp — fcp/(27r) ~ 3.092 x 10~^* Hz. The devices quoted in [10l[TP roughly probe regions between 60 and 110 
Hz while the ones of Ref. [T2|, IT3] are expected to work in the MHz region. 
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The idea pursued in this paper is, in short, the following. If the initial state of the 
gravitons is mixed, then it is also plausible that all the other fields (both conformally and 
non conformally coupled) will be in a mixed state characterized either by the same or by a 
similar comoving temperature. Among all the possible fields with different spins, particular 
importance should be attributed to the non- conformally coupled species. To address this 
broad problem it is useful to consider the simplest non-trivial situation where the initial 
state of the fluctuations is given as the direct product of the quantum state of the gravitons 
and of the other quantum states of the matter fields which will be assumed to consist of a 
collection of phonons (i.e. non- conformally coupled scalar fields): 



where the ellipses stand for other matter fields eventually present in the system. Both 
the gravitons and the phonons will be characterized by a density matrix whose statistical 
weights define the correlation properties of the initial state, as it will be specifically discussed 
hereunder. 

The secondary graviton spectra induced by a mixture of phonons are the main theme 
of the present analysis which may serve as a useful approximation to more realistic (but 
also necessarily more complicated) situations. The secondary graviton spectra, unlike their 
primary counterpart, are sensitive to the second-order correlation properties of the initial 
quantum mixture so that the semiclassical theory is argued to be generally inadequate. In 
section [2] the quantum treatment of the anisotropic stress will be outlined. In section [3] 
the discussion focuses on the relation between the secondary spectrum and the correlation 
properties of the initial state. In section H] the secondary spectrum is deduced in the realistic 
situation. Section [5] contains the concluding considerations and some perspectives for future 
developments in the light of forthcoming satellite experiments. 

2 Quantum anisotropic stresses 

A simple way to pin down the primary contribution to the graviton spectrum is to consider 
the limit where the gravitons are decoupled from the phonons but not from the initial mixed 
state. The secondary contribution depends instead on the contribution of the phonons to the 
total anisotropic stress which must be viewed, for the present purposes, as a field operator. 
Defining three mutually orthogonal directions as ki = ki/\k\, rhi = mi/\m\ and h = ni/\n\, 
the two polarizations of the gravitons in a conformally fiat background are: 




(1.2) 



4f\k) = {rhiThj - nifij), 



e\f\k) = {rhihj +nirhj). 



(2.1) 



From Eq. it follows that elj' e\j 



26 \x'. The sum over the polarization gives instead: 



J2 el^l{k) = Pmi{k)Pnj{k) + Pmj{k)Pni{k) - Pij{k)Pmn{k) 



(2.2) 



A 
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where Pij{k) = [Sij — kikj) is the traceless projector. The primary spectra of the rehc 
gravitons can be deduced from the following action: 

S=^l d'xV^^d^ hi d^h), (2.3) 

where g^j^ denotes a conformally flat background metric0 and £p = VSttG = 1 / Mp . In what 
follows we shall distinguish (especially in the last section) between Mp and Mp = 1.22 x 
lO^^GeV = a/Stt Mp. In the action (12. 3p hij denote the traceless and divergenceless modes 
of the geometry written in the form g^y^x.r) = g^i^{r) + Stg^u{x,T) where 6tgij = —a? hij 
and 5g'^^ = k"^ /a?. By taking the functional variation of the action (12. 3 p with respect to hi 
the equation of motion reads 

h'lj + 2m[j - V^h.j = 0, (2.4) 

where g^,^ = a^{T)r]^y. In Eq. (12. 4 p the prime denotes a derivation with respect to the 
conformal time coordinate r and T-L = (In a)' = a H where H is the Hubble rate. 



2.1 Primary spectra 



In terms of the two polarizations of Eq. (12. ip the mode expansion of the field operator 



hij {x, t) becomer^' 



h\]\x, t) 



(27r)3/2a(r) 



ik-x 



where and fk,\{,T) obey, respectively, 



f> 



kX 



0. 



(2.5) 



(2.6) 



The sum over A appearing in Eq. ( 12. 5p runs obviously over the two polarizations defined in 
Eq. ( 12. ip . The field operator of Eq. (12. 5p can also be represented, in Fourier space, as 

1 

- (n t\ = 



(27r)3/2 

V2ep 



£xhff{x, r) e''^-' 



(2.7) 



To compute the primary spectrum we need to calculate the expectation value 



hif{q^,r)h(^l^{q„T'm = TT 



h'^j\qur)h<;^J{q„T') 



(2. 



''Greek letters are used to denote four-dimensional indices; lowercase Latin characters denote spatial 
indices; uppercase Latin characters denote the different phonon indices. 

^The superscript in Eq. ()2.5p stands for the primary contribution which is the one obtainable by setting 
the phonon contribution to zero. 
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The first equality of Eq. fl2.8p follows from Eq. f ll.2p while the second equality prescribes 
that the expectation value must be computed by using the appropriate density matrix for the 
gravitons. We postpone to section [3] a more specific discussion of the density matrix itself. 
For the results of the present section the form of the density matrix is immaterial: what 
really matters is only the averaged multiplicity of the initial state and not its dispersion. 
The expectation value defined in Eq. fl2.8p can be written, in more explicit terms, using Eq. 

2il 



A A' 



2r^^ EEl /giA(^)/4A'(^04^^(^i)e^^2('?2)(agi,A42,A') 



+/;A(r)^AK04'V^i)eL'2(-^2)(aL,-.,,a_,-.,,v)}. (2.9) 

Using the Kroeneker deltas over the polarizations and the Dirac deltas over the momenta 
Eq. (12. 9 P becomes, for r — t- r' 

m^\q^,T) h^^i^iq^rm = ^ S^.^Qi) \U{r)? [2n(gi) + l]5^'\q, - q,). (2.10) 



a2(r) 



In Eq. (I2.10p the following auxihary sum Sijmn{k) = J2xe[j\k) el^l{k)/A has been defined; 
recalling Eq. (12. 7p . n{qi) denotes the averaged multiplicity of the initial state 



n{<li)^xx' -Q2) = (^l4i,Aa92,A'|V') = Tr 



Pgravitons ^(j\,A^?2>A' 



(2.11) 



In Eq. (I2.10p the mode functions of each polarization obey the same equation, i.e. where 
fqi® = = fq^. Recalling now that, from Eq. (12. 2p . Sijij{k) = 1, the primary power 

spectrum P^^ (g, r) arises from the first-order correlation between the graviton field operators 

m^\q,T)h§^\p,Tm = ^P^'\q,r)S^'\q-p), (2.12) 

4/2 3 

Pi'\Q,r) = -^^Am\'[2n{q) + 1]. (2.13) 



In the limit n(g) — )■ 0, using the results the mode functions calculated during a quasi-de 
Sitter stage (see, e.g. section H] Eqs. (I4.2p -( H731) ). Eq. (I2.13p gives exactly the properly 
normalized power spectrum of the tensor modes of the geometry of infiationary origin. 



2.2 Energy densities 



The energy-momentum tensor of the gravitons can be derived from the action (12. 3p by 
varying with respect to 'g°'^ . As argued by Ford and Parker [26] the resulting expression is a 
reasonable (but not unique) quantity to use for the calculation of the spectral energy density 
of the gravitons; the energy-momentum tensor is given, in this case, by: 



44 



1 

2^ 



(2.14) 
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The energy density of the gravitons is therefore obtained by averaging the (00) component 
of Eq. f l2.14p over the appropriate (mixed) state defined introduced in Eq. (11. 2p and aheady 
mentioned earher in this section, i.e. pcw = ('^I'To'lV')- Defining gq = and recalhng Eq. 
(12. lip the energy density of the gravitons is 

Pow = ^J d\nq£-^[\gq{T)\' + {q' + H')\f,{r)\' - H[f;{r)g,{T) + /,(r)^*(r)]} [2rl(g) + 1]. 

(2.15) 

By measuring Eq. (I2.15P in units of the critical energy density pcrit the critical fraction of 
energy density stored in the gravitons is obtained: 

1 dpGW 



PGW = (^I'T^IV'), 



crit 



r,2 f2 



(2.16) 



^Gw(g,^) 

Using Eq. fl2J5|) . Eq. ( l236|) becomes 

O , , Q'P^\Q.r) f^ , k(r)P + ^^|/,(r)p-?/[/;(r)^,(r) + /,(r)^*(r)] ^ 
^Gw{g,r)- 24^, |1 + -^^jj-^^ 1. (2.17) 

Instead of deriving the energy density from Eq. (I2.14p . it is equally plausible to use the 
energy-momentum pseudo-tensor [27] appropriately generalized to curved backgrounds [281 
[29] . By defining the energy- momentum pseudo-tensor from the second-order variation of the 
Einstein tensor, the expression for pcw is given by [301 l3T] 



1 



Pgw 



H{^\drhkih''\ij) + -md^mhMd"'h'''\^) + {^\drhMdrh''m 



(2.18) 



implying that ficw (<?>''") is now given by 

,2 p(l) 



o , , <l'P^\q,r) f_^ , l^7.(r)p-7?/^|/,(r)p + 3?^[/;(r)^,(r) + /,(r)^*(r)] 

""^^^^'"^ = -m^v + m^w 

(2.19) 

The difference between Eq. (12.171) and Eq. (I2.19P only appears in the correction to the 
leading results. By expanding Eqs. (I2.17P and (12.191) in the limit T-L/q the results are, 
respectively: 



^Gw(g,i") 



12V? 



12^2 



r 



1 



2g2 



(2.20) 



where we used thacl 5'g(r) ~ ±iqfq{T) for T-L/q <^ 1. In the opposite limit (i.e. q/V. <^ 1) 
the approximate solution of the mode function reads 



/,(r) ^ A(g)a(r) + B(g)a(r) 



dr' 
^2p)^ 



gq{T)c^nU + 



m 

a(r) ■ 



(2.21) 



^The following two statements are derived from the approximate solution of the equation for the mode 
function reported in Eq. (|2.6p . 
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For instance, during a quasi-de Sitter stage of expansion the term going as B{q)/a{T) is 
suppressed and the expressions of f^cwl^,''") coincide up to terms of order {q/l-lY as it can 
be checked by inserting Qqij) ~ l-Lfqij) into the general expressions of Eqs. fl2.17p and 



2.3 Anisotropic stresses 

The spectra obtained so far hold provided the induced anisotropic stress is strictly vanishing. 
In the opposite case Eq. (12 ■4p inherits a source term; by perturbing the Einstein equations 
and by focussing on the modes which are traceless and divergenceless we shall have 

1 



hi " + 2Hhl ' - V^hl 



(2.22) 



where 5tRl denotes the fluctuation of the Ricci tensor and where 11^ must be both diver- 
genceless and traceless, i.e. 

fi " + 2UU ' - V% = -2 a' fpfli , n: = dj fli, (2.23) 

where the classical fields have been promoted to the status of quantum operators (i.e. hij — )■ 
hij and 11] — ?■ 11^ ). The problem we ought to address is to obtain an explicit form of the 
anisotropic stress flj^ in Fourier space. By taking the Fourier transform of both sides of Eq. 
(12:231) we get 

+ 2Hhi, + q%j = -2lla\T)%{q, r), (2.24) 

where tiij{q, t) can be appropriately expanded in on the basis of the two polarizations of the 
graviton defined in Eq. (12. ip : 

n.,(g, t) = Y. ef{q) flA(g, r) = e® (g) fle(g, r) + e|(g)fl^(g, r), (2.25) 

A 

where 

r) = \Y. eZiq) vrL1](g, r), fl«(g, r) = 1 ^Zi^) ^i^M r). (2.26) 

^ A=l A=l 

The rank-to tensor vr^)(g, r) appearing in Eq. (I2.26P is in fact a convolution of two field 
operators: 

1 f 

TT^n ^) = (27r)3/2 a2(^) J (q ' P)0A{q - P, t) 0a{p, t) , 

C^^^{a,h) = arahn-'^-^ 5jnn, S = q - p, b = p. (2.27) 
In practice Trf^'' (g, r) is the Fourier transform of 

1 



-^ij {x,t) = 



(2.28) 



which is defined from the canonical energy-momentum tensor of the phonon field. In Eq. 
fl2.27p the field operators are defined as 



0Aik,T) 



-'kA 



kA 



e + Qa- 



kA 



0. (2.29) 



Eq. fl2.29p is derived by using a phonon action containing also the coupling to the Ricci 
scalar ESI, i.e. 



1 



AT 



^ = t; V-9 d'^x g°''^daipAdi3(pA -aRip\ 

^ •' A=l L 



(2.30) 



where R denotes the Ricci scalar. If a = 1/6 the coupling is conformal; if a 7^ 1/6 the 
coupling is non-conformal. We shall assume, in agreement with Eq. (12. 28 p . that the M 
scalars are minimally coupled (i.e. a = 0). Still it is useful to look at the general form of Eq. 
f l2.29p whose general solution, in the case of a quasi-de Sitter (and spatially fiat) background, 
is given in terms of Bessel functions. The dependence on a generic a can be easily included 
in the discussion and it amounts to a different Bessel index characterizing the evolution of 
the mode functions Fk^A- Along the same line, there is the possibility of including a mass 
term for the various phonons and this will result in a term m\o^ inside the squared brackets 
of Eq. (12. 29 p . During a quasi-de Sitter stage of expansion the mass contribution introduces 
a further modification of the Bessel index since m\a^ ~ (m^/if^)r~^. The fields (^a are 
fully inhomogeneous and do not have any background component. Their masses, if present, 
are assumed to be smaller than the infiationary curvature scale, as it will be discussed below. 
The set-up defined by Eq. (I2.30p is not exhaustive and it is only aimed at illustrating, in 
a simplified situation, how secondary graviton spectra can be produced and enhanced by 
stimulated emission. 

Using Eq. (I2.27P into Eq. (I2.26P we can get a more explicit form of n©(g, r) and il^{q-, r): 



n®(g,r) 



n®(g,r) 



1 



2(27r)3/2 a2(r) 



d^pCQ){q-p,p)0Aiq-p,r)0Aip,r) 

A 1 <^ 



1 



2(27r)3/2 a2(r) 



J2 J d^pC^iq-P,P)^A{q-p,T)<pA{p,T), (2.31) 



where, from Eq. (12. ip . 



C®(a, h) = {m ■ a){m ■ h) — {h ■ a){n ■ b), C^{a, h) = {fh ■ d){n ■ h) + {h ■ d){m ■ h). (2.32) 



Since, in Eq. (I2.3ip . a = q — p and b = p all the scalar products m ■ q and n ■ q vanish given 
that m, n and q are all orthogonal. As a consequence the final expressions can be simplified 
even further: 



C^{q-p,p) 



\q)PmPm C(^{q-p,p) 



\q) PmPn- 



(2.33) 
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Going now back to Eq. fl2.24p the evolution equation for each polarization can then be 
written as 



(2.34) 



We now have two complementary choices. Either we go on by representing the anisotropic 
stress on the basis of the tensor polarizations, or we go back to the tensor indices. In the 
second case, using Eq. f l2.34p we obtain: 



n.,(g,r) = eg)(g)fl,(g,r) = --^l-^J^J d'pSlf\q)p^ p-Mq-p,T)^Aip,T) 



(2.35) 



where the sum over the polarizations has been performed with the help of Eq. ( 12. 2p . i.e. 

S^MP) = P'"p"Eel?(^)4'i(9)=45.,„„(g)p'"p" 

A 



P 



{q-p)' 



5r 



1 + 



(g -p)^ 



m , ^PiPj 

2 ~'~ 2 



APi^j + qiPjJ^TY' 



q^p^ 



(2.36) 



It is simple to check that, indeed, Eq. f l2.35p is traceless and divergence-less fl* = gjll* = 0. 
To avoid confusions, note that Sij{q,p) = Ap^ p^ Sijmn{q) in agreement with the definition 
given right after Eq. (12.101) . 



2.4 Secondary Spectra 



Inserting Eq. f l2.35p into Eq. fl2.24p . the resulting expression can be simplified by introducing 
an appropriately rescaled anisotropic stress 



AT 



{A) 



h'l^ + 2Hh,, + q^K^ = 4 Q,, (g, r) , Q,, (g, r) = ^ Qf^ 

A=l 

Q\f{Q,^) = v,/o / d^pS'lf\q, p)0A{q-p,T)0A{p,r). 



(27r)3/2 

By rescaling hij, Eq. fl2.37p is further modified as 



(2.37) 
(2.38) 



fiij = ila{T) Qij (g, r) , /i^ (g, r) = a(r) hij (g, r) 



The solution of Eq. flOOD is 



i^iM r) = fifhq, r)+i'U a(0 G,{^, r) d^, 



(2.39) 



(2.40) 



where 

f (f) f*(r) - f*(f) f (t) 

r) = ' ^ > mo = m) gm - m) a.io = ^- (2.41) 

In terms of hij{q, r) the solution of Eq. (12 .40^ becomes: 

h^r) = h^\q,r)+il [G.^i^r) QiMOdi, G,{i,r) = ^G,{i,r). (2.42) 

Equation (I2.42p leads to the total spectrum of the gravitons in full analogy with the calcu- 
lation of the primary contribution reported in Eq. (I2.12p . i.e. 

{Mk,{q,T)hl{p,Tm = ^P^{q,r)6^'\q-p}, (2.43) 

where now Pt(?5 t) = Pt\q^ + Pt\(1j ''') denotes the total power spectrum written as the 
sum of the primary contribution obtained in Eq. f l2.12p and of the secondary contribution 
P^\q^ t) S^'^^q — p) which is explicitly given by the following expression: 

qHI. ^ ^ 



27r2 



E E r^^i r di^G,{i,, r)G,i^2,r) (v&l sSf (g, 6) SSf (-P, 6)1*)- (2-44) 

A=l B=l '^^ 



Note that in Eq. (I2.44p we used that Q|j(p, = Qiji^P^O- must be also appreciated that, 
as a consequence of the normalization of the Wronskian in Eq. (I2.4ip . we have Gg(^,r) = 

g;{^,t). 

To conclude this section it is appropriate to mention that the relation holding between 
the primary spectrum and the spectral energy density f2Gw(?) t) can be obtained also in the 
situation where the secondary spectrum is present. In particular it can be shown that 



(2.45) 



where the power spectra are evaluated for typical scales larger than the Hubble radius (i.e. 
qr -C 1) prior to matter radiation equality. The function f{q, ggq) parametrizes the transfer 
function of the tensor spectra across matter-radiation equality [9J 

/(g, geq) = 1 + ci (^^) + h (^) ci = 0.5238, h = 0.3537, (2.46) 

where 

.e, = ^ = 1.288 X 10-^ ('^] ( \ Hz; (2.47) 

^ 27r V0.1334/ V4.15 X lO-sy ' ^ ' 

VLyio and ^ro denote, respectively, the present critical fractions of non-relativistic matter and 
radiation. 
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3 Quantum correlations 



The secondary power spectrum of Eq. fl2.44p is determined by the expectation value 

6^(gi,g2;6,6) = i^lQ^fiqu^i) QSf (-^2,6)1*), (3.1) 

whose exphcit form is given by the following expression: 

^(gi,92;^l,6) =/ (^271)^'^ I (27r)3/2 ^iA)(^^'^P^) S'^B)i?2j2)MAB{qi, Pi] 52, ^2)- (3.2) 

The explicit expression of A^a_b('?i, Vi] <l2, P2) appearing in Eq. (13. 2 p is quartic in the phonon 
fields: 

MAsiqu Pi, <?2, P2) = (*|<^A(gi -Pu ^i)<^a(pi,^i)<^U^2 -P2,6)^Up2,6)I*), (3.3) 

and (pc{k,T) = h^^FkciT) + 6^g^-F^*(^(r). In Eqs. (13. ip and (13. 3 p the average over the 
quantum state determines the properties of the final spectrum. 

The primary graviton spectrum p!^\q^T) depends on the averaged multiplicity of the 

(2^ 

initial state. The secondary graviton spectrum Pr^ {q,^) depends on the expectation value 
of two anisotropic stresses i.e. four field operators. In terms of the statistical properties of 
the initial state this means that secondary graviton spectra are not only sensitive to the av- 
erage multiplicity of the initial distribution but also to its variance. There is here an analogy 
with quantum optics [22] where first-order correlation effects (i.e. Young interferometry) 
are unable to define completely the statistical properties of the light source. Second-order 
correlation effects are instead probed with Hanbury-Brown-Twiss interferometry where, in- 
stead of correlating the field amplitudes the correlation is between field intensities which 
are themselves quadratic in the field amplitudes. Similarly, primary and secondary graviton 
spectra depend, respectively, on the average of field amplitudes and field intensities. In this 
connection, two issues need to be addressed: the problem of the ordering of the operators 
and the statistical properties of the initial quantum state. 

3.1 Ordering of the operators 

Should we impose a specific ordering to the operators appearing in Eq. (13. 3p ? If we do 
not impose any ordering, the averages of the creation and annihilation operators appear in 
various combinations. Following the standard habit of quantum optics [32] and suppressing 
all the indices it is useful to consider the following normalized expectation values in a simple 
quantum mechanical framework where only one mode of the field is considered: 

jm_ Cb^brbb) j^,_ {mh^) j,,_ {hhrwh) 
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where [6, 6^] = 1. The superscripts of Eq. fl3.4p denote, respectively, the normal, antinormal 
and symmetric ordering. The averages are taken over a mixed state characterized by a 
density matrix 

oo oo 

P=Y.'^n\n){nl E^n = l- (3-5) 

n=0 n=0 

Using the commutation relations Eq. (13. 4p can be written without mentioning any specific 
property of Vn as 

9^^' = g'""' + A + 9^'' = a'""' + 7^v a'""' = i + (3-6) 

{n) (n)^ (n) (n)^ 

where 

(n) =TT[pb^b], (n^) = TT[pb^bb^b], = {n") - {h)\ (3.7) 

Consider now two examples associated with the density matrix of Eq. (13. 5p : a Poisson 
mixture (denoted by V^^) and a Bose-Einstein mixture (denoted by V^^^)] their explicit 
expressions are given, respectively, by: 

^r = ^e-, ^r^ = ^(^)'^- (3.8) 

n\ n + l\n + 1/ 

In the case of the Poisson mixture g'^^'^ = 1 since, in this case, (n) =n and = n as it can 
be expected from the well known properties of the Poisson distribution. The antinormal and 
symmetric correlation functions g^^^ and g^^^ coincide, approximately, in the large n limit 
with the normal-ordered expression g^^\ 

In the case of the Bose-Einstein mixture g^^'' = 2 and, as in the Poisson case, g^^^ ~ g^^'^ 
for n ^ 1 since, from P^^^^ of Eq. (13. 8p . = rf + n. Again the latter result is implied by 
known properties of the Bose-Einstein (or geometric) distribution. Let us finally consider 
the interesting example of the negative binomial mixture where the weight is given by: 

r(C)r(n+l) \n + (J \n + CJ ' ^ ' ^ 

Equation (13. 9 p corresponds to the distribution followed by ( Bose-Einstein distributions all 
emitting at the same temperature. In the case of Eq. (13.91) 

^(^^ = 1 + 7, ^ = 1 + 7, (3.10) 

with g^^^ and g^^^ coinciding in the large n limit. The examples of Eqs. (I3.8p - (l3.10p sug- 
gest two complementary conclusions. The first conclusion is that since the normal-ordered 
expectation values dominate the normalized correlation functions, it is useful to employ the 
normal ordering as long as the interesting physical limit is the one for n ^ 1. The second 
conclusion is that the semiclassical approximation fails in the case of a Bose-Einstein initial 
state and, more generally, whenever — n^O. In quantum optics it is natural to impose 
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the normal ordering in the higher-order correlators since the detection of light quanta (i.e. in 
the optical range of frequencies) occurs by measuring a photo-current, i.e. a current induced 
by the absorption of a photon [32]. In the present case by normal-ordering the expectation 
values of two intensities (i.e. four field operators) we implicitly demand that when the av- 
eraged multiplicity of the phonons goes to zero the secondary contribution to the graviton 
spectrum vanishes. This is consistent with the requirement that, ultimately, the interesting 
physical regime involves typical wavelengths much larger than the thermal wavelength where 
the averaged multiplicity of the (thermal) phonons is large. 



3.2 The density operator and the expectation values 

In the quantum mechanical examples discussed so far the dependence on the momenta has 
been neglected but their presence is essential. The expectation values are evaluated using 
the appropriate density operator. Recalling Eq. (11. 2p the density operator can be written 
as p = Pgravitons ® Pphonons- In the Fock basis the gravitons and the photons will then be 
characterized, in principle, by different statistical mixtures since the two density matrices 
span two different Hilbert spaces. The notations established in section |2] (see, e.g. Eqs. (12. 7p 
and (12.291) ) imply that the gravitons are characterized by creation and annihilation operators 
conventionally denoted by at^ and ag^ while, for the phonons, the notation bl^ and has 
been employed. The density matrix for the gravitons is therefore given by 

Pgravitons = E E ^{^r""^^^ = 1, (3.11) 

{m} {m} 

where = 1^^^^) ll"^fc2A) W^itsx)--- ellipses stand for all occupied modes of the 

field. Similarly the density matrix for the phonons can be written as: 

Pphonons = E ^S"^ E ^S"^ = 1- (3-12) 

{n} {n} 

The P{n} (either for the phonons or for the gravitons) will be taken in the form0: 

p,„, . n .:':^;"^' f^)'" (^T\ (3.13) 



, r(CE)r(nj + 1) V5Tj + (fj Vnj + c, 

which is the field theoretical generalization of the quantum mechanical example of Eq. (13. 9p . 
In the limit (j^ ^ 1 the distribution of Eq. (I3.13P becomes a Bose-Einstein distribution: 

— "fe 

n..i^p{ff'=n (i/i (3.14) 



^To avoid the pedantic repetition of the same formula with shghtly different subscripts the indices denoting 
the polarizations of the gravitons (i.e. A) and the different phonon fields (i.e. A) have been dropped. 
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Using Eqs. fl3.1ip and fl3.13p the expectation values are computed with the help of the 
density operator: 

(3.15) 



Pphonons ^q^p^q' ^p' 



Using a stenographic notatioii^ where a each letter corresponds to a momentum the previous 
expectation value can be written as 



(6| b] bk bi) = {b\ b] bi bi)6ij djk^ik + 



{b\b] bibj) 5ik5je[l - %] + {bjbjbj h) 5it5jk[^ - 5i 



(3.16) 



In the term (S| bj hi hi) there are four different operators all acting on the same momentum; 
in the second class of terms (i.e. (6| hi hj) and (hi 6] hj hi)) the momenta are paired two by 
two (as the corresponding deltas indicate). Using Eq. f l3.13p we shall have that 



{hi hi hi hi) = mi{mi - l)Pj(mi), 

mi 

{hi h] hj hi) = Pi j ' "^i ) ' 



blh] 



bj bi bj) = Y Pi ji^i^ ) 

rrii, raj 



(3.17) 



where Piirrii) = Pk^mk,), Pij{mi,mj) = Pk^mk,) Pk^irUk^) ■ The relations of Eq. (I3.17p 
follow from the definition of P{m}) i-e. P{m} = Pkiif^^ki) Pk2{f^k2) Pk^ifT^ki) ■■■■ where the 
ellipses stand for the product over the various momenta. Inserting Eq. fl3.17p into Eq. 
f l3.16p we get, after simple algebra, 



{hi h] hk hi) = nl(^ - l^5ij 6jk 6ke + rii rij 



^ik^je + ^ii^jk 



(3.18) 



In the limit Q ^ 1 for every mode of the field the contributions of the terms of the type 
(hi hi hi hi) exactly cancels; the standard rules of evaluating correlators in thermal field theory 
is quickly recovered [33]. If, however, Q I the first contribution in Eq. (13.180 does not 
vanish. The same procedure leading to Eq. ( I3.18P can be used to compute, with the due 
differences, all the other expectation values arising when the ordering is either anti-normal 
or symmetric. Just to give some examples consider, for instance, the following results: 



(hi hj hlh\) = n^i^— -l) 5ij 5jk Sk£+ (ni + l) {rij + 1) 



{hh]hlh,) 



'^i ( ^ - 1 ]^ij^jk5kt + rii>{ni + 1) 



^ik^je + ^ie^jk 
SijSke + SikSej 



(3.19) 
(3.20) 



^The operators bi and &| denote the annihilation and creation operators related two generic momenta, i.e. 
for instance and bt ; furthermore, following the same shorthand notation, Sij denotes the delta functions 
over the three-momenta (i.e. S^^p). 
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3.3 Secondary spectrum 

Using the notations ki = qi—pi, and k2 = 52—^2 the expectation value of Eq. (13. 3p becomes 
MABiquPi; 92, P2) = {^\blr.bl^^,Ab^k,BK^2B\^)F:^A^^)F;^A^^^ 



H'^lbi J, &k B K A hp,A\^)Fl^{^,)F;^s{^2)Fk,A{^i)F,,A{^i) 

H%K^,A\sKAK^.B\^)Flj,{i,)F,,B{i2)F,,A{il)^^^^ 
+ ('^\^U.A\B^^V^^^^-P^BmFl^s{i,)F^ 



(3.21) 



Recalhng that [%^, ^^5] = S^'^\k — p) 6ab, the resulting expression for the secondary power 
spectrum becomes then: 



g2p2 

J^2fi4 ^7-2 



1 2 



- 1 



647^5 



^n/^ (f " ly^' /Jrf6G,(ei,r) G,(6,r)W(g;6,6), (3-22) 



where V denotes the fiducial normalization volume and the other functions are defined as 



J^{k; ei, 6) = Fk{^^) F*{^2) + F*i^,) F,(6), 

= [F'{^i)f;\^2) + F*\^,)F'{^2)]. 



(3.23) 
(3.24) 



The second term in Eq. (I3.22p vanishes in the Bose-Einstein limit (i.e. (g — )■ 1) but it 
vanishes also in the infinite volume limit. This term is per se interesting but its analysis is 
not central to the present discussion. So, in what follows, we shall just assume that Cq = I 
which means that the phonons will be characterized by a single Bose-Einstein mixture. 



4 Bose-Einstein enhancement 

To assess the Bose-Einstein enhancement of secondary graviton spectrum we must go back 
to Eqs. (12.411) and (I2.42p and compute explicitly Gfc(^, r) and Gk{C,, t). For the estimates 
of the integrals appearing in the expressions of the secondary graviton spectra it will also 
be necessary to deduce more explicit forms of Eq. (13.231) . Equations (12.61) and (I2.29P can be 
solved using standard techniques [Ml [351 ES] in the case of quasi-de Sitter space-times. 
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4.1 Explicit expression for the Green's functions 

From Eq. fl2.4ip the explicit expression of Gq{^, r) during a slow-roll phase is given by 



TT 



Yu{-kO M-kr) - Y,{-kr) M-k^ 



2(1 -e)' 



(4.1) 



where e = —H/H"^ is the slow-roll parameter. Equation (14.1 p are derived by first solving Eq. 
(12. 6 p and the related equations for the tensor mode function^ 



9i{k,T) 



N 



'2k 

drfk 



k 



N 



(2) , 



-x) + 



1 -2z/l 



-X] 



Hi'\-kr) 



(4.2) 
(4.3) 



where H^\z) is the Hankel functions of first kind [M]. Notice that Eqs. (14. 2 p and (14. 3 p imply 
a constant Wronskian, i.e. W{S,) = i- For power-law inflation we have a(r) = (— T*/r)^^ 
and u = {P + 1/2). In the latter case the rescaled Green functions are expressible as 



2k 



T 



Y,{-ki) -M-kr) - Y^i-kr) J,{-ki) 



In the limit e ^ 1 we have that Gk{C,, t) is, approximately, 

1 



Gfc(e,r) 



k{^ - r) cos m - r)] - {k'^T + 1) sin [k{^ - r)] . 



(4.4) 



(4.5) 



It is also possible to deduce a more general expression for the Green's function valid across 
the transition between inflation and radiation. The expression for the mode function is 

1 



fk{r) 



'2k 



-r - rj I 1 - — |e 

KT 



kT^ 



1 



2Pr2 



-ikr 



1 



2Fr2 



(4.6) 



where r* denotes the transition time between the inflationary phase and the radiation- 
dominated phase and 6{z) denotes the Heaviside theta function with argument z. To derive 
Eq. (14. 6 p it is necessary to use an expression for the scale factor which is continuous across 
— such as 



a(r) 



-r* — r 



r + 2n 



(4.7) 



When the mode function is given by Eq. (14. 6p . the rescaled Green's function is 



G,i^, r) 



-T - n 



^ g(e - r) cos [g(e - r)] - {q^r + 1) sin [q{i - r)] 



r + T, 



(4.8) 



^The normalization N appearing in Eqs. (|4.2p and (j4.3p should not be confused with J\f counting the 
phonon fields. 
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To estimate the secondary graviton spectrum the functions defined in Eq. fl3.23p must be 
computed; the function F{k] ^i, ^2) is given, in exphcit terms, by 



2a(6)a(6) 



(4.9) 



where the Bessel index /i depends on the couphng of the phonons, i.e. 7^ 1/2 if the couphng 
is not conformal; fi can also contain an exphcit dependence on the mass of the phonon as 
aheady remarked in connection with Eq. (12.29^ . For sake of simphcity the mass of the 
various phonon species will be assumed to be always smaller than the de Sitter curvature 
scale. In the limit k^i ^ 1 and k^2 ^ 1 the function J^{k; ^1, ^2) can be expressed as: 



^'^''^'^ A;7ra(6)a(6) V 4 J a(6)a(6) ^'^r I 4 



3/2-^1 



(4.10) 



Instead of performing the integrations over 6 and 6 iii Eq. fl3.22p it is useful to introduce 
the following set of rescaled variables 



p k 
q q 



xi=q^i, x2 = qC.2, y = qr. 



(4.11) 



For instance, with the rescalings of Eq. (14. lip the explicit expressions for the functions 
Gq{^,T) and J^{p; ^1,^2) are, for r < -r*: 

(xi - y) cos (xi -y)- {xiy + 1) sin (xi - y) 



Gq{ii,T) = Gg{xi,y) 

Gqi^2,r) = Gq{X2,y) 



qxi 



(x2 - y) cos (x2 - y) - {x2y + 1) sin (x2 - y) 
qxi 



J^iPi, ^1, 6) 

Hki; 6, 6) 



(a^xiX2 + 1) cos [a{xi — X2)] — a{x2 — xi) sin [a{xi — X2)] 

p a{^i) a{^2) ct'^ X1X2 
(/3^xiX2 + 1) cos [l3{xi - X2)] - I3{x2 - xi) sin [/3(xi - X2)] 



ka{^i) a(6) 13'^ xi X2 
where we remind that, according to the notations of Eq. (13.211) . k = \q—p\. 



(4.12) 



4.2 Explicit expressions for the secondary spectrum 

With the explicit formulas derived so far, the secondary spectrum can be expressed as: 

P?\q^^) = ^Jd'P Jl^^ d^2 C{q, p- 6, 6, r) 7^(g, p), (4.13) 

where the functions £(g, p; 6; t) and 7^(g, p) are: 

>C(g, p; 6, 6, r) = G,(6,r)G,(6, T)^(|g-p1;ei, 6)-^(p;ei, 6), (4.14) 



T^iQ^p) = Q P n{\q—p\)n{p) 



g2p2 



(4.15) 
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Let us now consider, as an example, the case of minimal coupling of the phonons leading to 
/X ~ 3/2. In this case the integrals over and ^2 can be performed exactly by breaking the 
integration for r < — r^, and for r > — r*: 



r d^i r ^6 Ag,P;ei,6,r) = ^^^^Ji4^, (4.16) 



where the function Z{t, tq, r*) is 

_ qro cos [g(r + Tq)] - sin [g(r + tq)] sin [g(r + rj] - gr cos [g(r + rj] 

-^IT, To, T*) — h 



qro qn 
qn COS [g(r + r*)] + sin [g(r + r*)] - g(r + 2r=^) 
q{r + 2n) ■ 



(4.17) 



The amplitude of each intensity (quadratic in the field operators) has been denoted by 
A = A/ip] from the latter relation it also follows that A = £pA is the dimensionless 
amplitude given in Planck units. The correct physical limit where the function Z{t, tq, r*) 
should be evaluated is the one where gr <^ 1, gr* ^ 1 and |ro| ^ |r*|. Recall, indeed, 
that Ti = —Tq corresponds to the initial integration time and — marks, in the sudden 
approximation, to the end of the quasi-de Sitter phase. After this limit the dependence on 
r drops, to leading order. Bearing then in mind the expression of TZ{q,p) of Eq. f l4.15p and 
putting all together we have that the expression of the secondary spectrum can be written 
as: 

P?\l)-^fp'''p'^i^\'^-A'- (4.18) 

167r^ J q\q — p\'^ i q p J 

The dependence on the conformal time coordinate has been dropped since we are considering 
here the secondary contribution for the modes which are larger than the Hubble radius prior 
to matter-radiation equality. This is the quantity which must be directly compared with the 
standard form of the primary spectrum (see section [5]). 

In Eq. f l4.18p A denotes the common amplitude of the phonon fields; A depends on 
the Hubble rate in Planck units during the quasi-de Sitter stage of expansion and will be 
specified later for notational convenience. Without fine-tuning A ~ A-r where An is the 
amplitude of curvature perturbations (see also section [5]). The most interesting limit for the 
present analysis is the one where n{p) ^ 1 and n{\q — p\) ^ 1. Thus n{p) and n{\q — p\) will 
be approximated as 

1 (Jt 

^(g) = e<?/^T - 1 - 7' (^-^^^ 

where q^ = T is the (comoving) thermal wavenumber; we shall also assume that n{q) = 
for q > q^. Using Eq. fl4.19p inside Eq. fl4.18p the secondary spectrum can be written as 

4^^(g) = ^^X(g,gT,A:p), (4.20) 
167r^ q 
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where kp = 0.002 Mpc ^ denotes the pivot scale aheady mentioned in the introduction and 
corresponding to a pivot frequency Up = kp/ (27r) = 3.092 x 10~^® Hz. The integral appearing 
in Eq. (OOj) is: 



I{q,qT,kp) 



(q-p)^' 



1 2 



q2p2 



\q-p\^ 

By making explicit the angular and radial integrations and we shall have that 

ix' - 1? 



X(g, gT, kp) 



2tt 



dx 



p dp 



(g2 + p2 _ 2qpxy ' 



(4.21) 



(4.22) 



where we used that \q — p\ = \/ q^ ^ p)^ — Iqpx. By explicitly performing the angular inte- 
grations in Eq. (14.221) . the explicit expression of the secondary spectrum becomes 



4'^ (9) 



A^M ( q 
1927r4 Vgx 



IT dp (A [Sp'^ - 2g2p2 _^ 3g 
P P I 



q3 p 



41 (2 



q^f (P + T 



p2 q4 



In 



\p-q\ 



\p + q\ 



(4.23) 

To estimate the last angular integral appearing in Eq. fl4.23p we divide the integration 
interval in two regions, i.e. kp < p < q and q < p < qT- In the first region we can expand 
the integrand in powers of \q/p\ < 1 and perform the integration over p. Similarly, in the 
second region the integrand can be expanded in powers of \p/q\ < 1. The final result for the 
secondary spectrum for kp < q < q^ can then be expressed as: 



52AfA^ f q 



315 TT^ Vgx 



1 



7_ 

26 V g 



U \ 3 

tvp 



26 Ut 



+ 



1 / g y 



13 VgT/ 



(4.24) 



Finally, as already remarked, for q ^ qT the secondary spectrum is exponentially suppressed. 



5 Concluding considerations 

Physical considerations suggest that the gravitons and the phonons have the same tem- 
perature. However different situations cannot be ruled out. In what follows the primary 
and secondary contributions will be compared in different limits. By focussing on the pri- 
mary contribution, let us insert Eq. (14. 2 p and Eq. (I4.19P into Eq. (I2.13p : in the limit 
|gr| = \q/{aH)\ ^ 1 the result is: 

^ 4^^"5r^(^)(l - e)-.(^)^-^oth (JL). (5.) 

which is nothing but the amplitude of the primary power spectrum for typical length-scales 
larger than the Hubble radius prior to matter-radiation equalitj0. As already mentioned in 
Eq. (14. ip . z/ = (3 — e)/[2(l — e)]. There are now different (but equivalent) ways of expressing 

^'^ In Eq. ()5.ip . the term (1 — e)^"^^ arises by eliminating r is favor of {aH)~^ since, within the slow-roh 
approximation aH ~ — 1/[t(1 — e)]. 
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the result of Eq. (15. ip . Recalling that ip = Mp^ the spectrum (15. ip at horizon crossing (i.e. 
q ~ Ha) becomes 



P¥'iq) = ^ =^ coth , (5.2) 



where the slow-roll relation 3MpH^ ^ V has been used and V denotes the inflaton potential. 
It is sometimes useful to express the normalization of Eq. (15. 2 p in terms of Mp = a/Stt Mp 
with the result that the prefactor [2/(37r^)](y/Mp) appearing in Eq. (15. 2p is modified and 
becomes (128/3)(V"/Mp). From Eqs. (15. ip and (15.30 . to first-order in the slow-roll approxi- 
mation, the standard tensor spectral index nx is 

nT = 3-2z/~-2e + C(e2), tt = 16e = -Snx, (5.3) 

where rx is the ratio between the tensor and the scalar amplitude in the concordance model. 
The relation between the tensor amplitude and the infiaton potential is not central to the 
present consideration so that the primary graviton spectrum can be usefully parametrized 
as ^ 

pW(g) = ^xff coth f-^y = r^An, (5.4) 



.kpj \2q' 

where = ^tAu is the tensor amplitude expressed in terms of the amplitude of curvature 
perturbations An at the pivot scale fcp = 0.002 Mpc~^. The values of the cosmological 
parameters determined using the WMAP 7yr data alone [H [2] in the light of the vanilla 
ACDM scenario ar 

(^]b, ^^c, f^de, /io, ere) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088). (5.5) 

If a tensor component is allowed in the analysis of the WMAP 7yr data alone the relevant 
cosmological parameters are determined to be [D E] 

(fib, f^c^^de, /io, ere) = (0.0430, 0.200, 0.757, 0.735, 0.982, 0.091). (5.6) 

In the case of Eq. (15. 5 p the amplitude of the scalar modes is An = (2.43 ± 0.11) x 10~^ 
while in the case of Eq. (15. 6 p the corresponding values of A-r. and of rx are given by 

y^Te = (2.28 ±0.15) X 10"^ rx < 0.36, (5.7) 

to 95 % confidence level. The qualitative features of the effects discussed here do not change 
if, for instance, one would endorse the parameters drawn from the minimal tensor extension 
of the ACDM paradigm and compared not to the WMAP 7yr data release but rather with 
the WMAP 3yr data release, implying, for instance. An = ^.itl'^ x 10~^, Us = 0.984 and 
rx < 0.65 (95 % confidence level). 



^^Following the standard notations r^c, and f^do are, respectively, the critical fractions of cold dark 
matter, baryons and dark energy; ere denotes the optical depth at reionization and rig is the spectral index 
of curvature perturbations; ho is the Hubble rate in units of 100 km/ [sec x Mpc]. 
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Let us now compare the primary and secondary contributions to the graviton spectrum 
in two different and extreme cases. In the first case the gravitons are not produced from 
stimulated emission but rather from spontaneous emission. This hmit is recovered from Eq. 
f l2.13p by setting n{q) — )■ 0. Conversely the phonons are produced by stimulated emission. In 
the latter case the primary spectrum, unlike the secondary contribution, will not experience 
any Bose-Einstein enhancement. Consequently, the total power spectrum will be given by: 

P,fe).^,ar%5!^m-'. (5.8) 



kpj 315 vr^ Vgx 

where only the leading term of the secondary contribution has been kept (compare with 
Eq. (14.241) ). In Eq. (15. 8 p the amplitude of the phonon field has been estimated by setting 
A ~ Ati. This choice is fully justified since this is what happens in the case of a quasi-de 
Sitter phase. 

To evaluate explicitly the figures in Eq. (15. 8p we still need an estimate of the ratio q/qT- 
Let us therefore proceed along this direction and remind the definition of A^max, i-e. the 
maximal number of efolds which are today accessible by our observations: 



ax 



i2..A.y^'i^..y^'{!^y''\!^)--"\ (5,9) 



where Hq = 100 /iq Mpc~^ km/sec is the present value of the Hubble rate and h^^l^Q = 
4.15 X 10~^. In Eq. (15. 9p we included the possibility of a delayed reheating terminating at 
a putative scale Hr possibly much smaller than the Hubble rate during inflation denoted 
by H; 7 controls the expansion rate during this intermediate phase. Note that Hr can be, 
at most, 10~'^^Mp corresponding to a reheating scale occurring just prior to the formation 
of the light nuclei. If 7 — 1/2 > (as it happens if 7 = 2/3 when the post- inflationary 
background is dominated by dust) A^max diminishes in comparison with the case when H = 
Hr- Conversely if 7 — 1/2 < (as it happens in 7 = 1/3 when the post-inflationary 
background is dominated by stiff sources) A^max increases. If Hr = H (or if 7 = 1/2) there 
is a sudden transition between the inflationary and the post-inflationary regimes and, in 
this case, we have approximately A^max — 64 -|- 0.25 In e. The value of A^max has then a 
large theoretical error, i.e., based on the previous considerations, 65 ± 15. In practice A^max 
is determined by fltting the redshifted inflationary event horizon inside the present Hubble 
radius Hq^. If the total number of inflationary efolds A^t is larger than A^max (i-e. A't > A^max), 
then the redshifted value of the inflationary event horizon is larger than the present value of 
the Hubble radius. The value of the thermal wavenumber is therefore given by 

A^max) 




where A/th denotes the total number of spin degrees of freedom in thermal equilibrium at the 
onset of inflation. From Eq. (I5.10p it is possible to argue for which values of the parameters 
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q < so that the possible enhancement of the power spectrum in Eq. (15. 8p is effectively 
realized. In this respect from Eq. (I5.10p it follows that 

' - O(10)e-("0f_^)<^f_^)''V (5.11) 



qt ' ' V/toMpc"V V2.2OXIO-V VO.l 

Choosing A^t ~ A^max and focusing the attention on the pivot scale (i.e. g ~ A;p in Eq. (15. Sp ). 
the total spectrum of Eq. (15. 8p becomes: 



PT(fcp)^^T+ :.vr . (5.12) 

315 TT* \qT/ 

Let us now suppose, in agreement with Eq. (15. lip , that k^/qi — 0{1Q~^) (at most) and 
M ^ (9(10^) (at most). In the latter case the first term in Eq. (15.120 must be smaller 
than 10^^*^ (according to Eq. (15.70 ) but the second term in Eq. (15.120 can be as large as 
C(10^^'^). Let us finally suppose that both the gravitons and the phonons are produced from 
stimulated emission. In this case Eq. (I5.12p is modified as 



P^{q) ~ At 



315 At Vgx 



(5.13) 



The primary tensor spectrum is also enhanced so that the secondary contribution is always 
smaller except in the limit of a very large number of phonon fields such as A/" > (9(10^). 
Note that Eq. (I5.13P holds for q < q^; in the opposite limit (i.e. q > gx) the secondary 
contribution vanishes while PT^q) — AT{q/kp)"'^ . If the total number of efolds is Nt ~ iVmax 
it seems plausible that the primary and secondary contributions could be directly bounded 
by forthcoming satellite missions such as Planck explorer. The phenomenological signature 
will be a power spectrum with primary slope (wt — 1) and with secondary slope going as —2. 

The same hierarchies between the primary and secondary spectra are reflected, according 
to Eqs. (I2.45p -f l2.47p . in the spectral energy density Qc-^i^q, r) whose explicit expression can 
be written, function of the frequency, as 

'*-(..r„) = [,.,(^)"%oth(£-)+5i^(-^)"]. (5.14) 

-,-15/^ ^O^RO \f 



= 3.94x10-1^ J • (5.15) 

V4.15 X 10-V V2.28 X 10-V ^ ' 



For u > z/max — GHz ilG\\/{u,To) is exponentially suppresseco In spite of possible observa- 
tional implications (which are notoriously difficult to assess in the case of the tensor modes 
of the geometry) the theoretical lesson we can draw from the present analysis is that the 
correlation properties of the initial quantum state can modify the tensor power spectrum 



^^For the fiducial set of parameters adopted above we have that i^max = (e/0.1)^/'*[y^7z/(2.28 x 
10~^)^/*][/iQriRo/(4.15x 10^^]"'^/'' GHz. These frequencies wiU be however not touched by the present cosider- 
ations. 
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at large-scales. This situation bears some analogy with the Glauber theory of optical co- 
herence: first-order coherence of light is controlled by the analog of the primary graviton 
spectrum while second-order coherence is described by the analog of what has been called 
secondary power spectrum. Consequently, the secondary graviton spectra, in spite of their 
detectability prospects, reflect the statistical properties of pre-inflationary mixed states and 
of their second-order correlations. 
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